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stabilized quantum gravity is equivalent to Langevin evolution (i.e. stochastic quan- 
tization) between fixed non-singular, but otherwise arbitrary, initial and final states. 
The simple restriction to a fixed final state at £5 — > 00 is sufficient to stabilize the 
theory. This equivalence fixes the integration measure, and suggests a particular 
operator-ordering, for the fifth-time action of quantum gravity. Results of a numeri- 
cal simulation of stabilized, latticized Einstein-Cartan theory on some small lattices 
are reported. In the range of cosmological constant A investigated, it is found that: 
1) the system is always in the broken phase < det(e) >^ 0; and 2) the negative 
free energy is large, possibly singular, in the vincinity of A = 0. The second finding 
may be relevant to the cosmological constant problem. 
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1 Introduction 



The "fifth-time" action is a general procedure, first proposed in ref. |]J, for sta- 
bilizing theories whose Euclidean actions are unbounded from below. There are a 
several such "bottomless action" theories which are of interest to physics; examples 
include D=0 matrix models, bosonic string field theory, and, especially, Einstein- 
Hilbert gravity. The fifth-time method has been applied by many authors to 
D=0 matrix models. In this article I will continue the study of fifth-time stabilized 
quantum gravity, begun in ref. ||. 

Euclidean quantum gravity is of interest for at least two reasons. First, there 
are the intriguing Baum-Hawking-Coleman arguments [|J for the vanishing of the 
cosmo logical constant, which are formulated in the Euclidean approach. Unfortu- 
nately, these arguments involve the semiclassical evaluation of a badly divergent 
path- integral, and are therefore rather suspect ||. Secondly, Monte Carlo simula- 
tions of quantum gravity cannot avoid the Euclidean formulation. In ref. |5] it was 
shown how the fifth-time prescription generates a new diffeomorphism-invariant ac- 
tion for quantum gravity, bounded from below, which leads to the usual Einstein 
equations of motion in the classical limit. Moreover, the stabilized Einstein-Cartan 
theory appears to be reflection-positive, at least in lattice formulation, and should 
be a good starting point for non-perturbative investigations. 

At the perturbative level, the stabilized theory flips the sign of the "wrong-sign" 
conformal mode at zeroth order, and is non-local at higher orders. This behavior is 
in agreement with two other approaches which aim at deriving the Euclidean theory 
from Minkowski space gravity. The first, due to Schleich (see also 0), solves the 
Hamiltonian constraints in the Minkowski theory prior to Wick rotation, and then 
reinserts redundant degrees of freedom to recover diffeomorphism invariance. The 
second approach, due to Biran et. al. || and Mazur and Mottola ||, argues that the 
kinetic term of the conformal factor should be Wick rotated like a potential term, 
when certain Jacobian factors in the integration measure are taken into account. 
In both bounded, non-local Euclidean action is generated. These studies 

show that the Euclidean continuation of Einstein-Hilbert gravity is not necessarily 
obtained by the naive replacement of the Minkowski action by the corresponding 
(unbounded) Euclidean expression. The advantage of the fifth-time approach is that 
the 5-dimensional action is local, relatively simple, and amenable to numerical and 
analytical methods, whereas the non-local stabilized 4-dimensional action, obtained 
in any of these approaches, is not known in closed form. 

A common objection to the fifth-time prescription is that it is rather ad hoc, at 
least as presented in the original treatment in [l]. The method generates a stabilized 
action which has the same classical limit and formal perturbative expansion as the 
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unstabilized theory, but there could be other methods to achieve those ends. In 
fact, for D=0 matrix models alternate stabilizations have been proposed [TI], Oj]. A 
better understanding of the physics underlying the fifth-time action would therefore 
be helpful in judging its merits. In the next section it will be shown that the 
fifth-time prescription is equivalent to stochastic quantization, with the constraint 
that the Langevin evolution not only begin at a fixed (arbitrary) initial state (at 
tr, — > — oo), but also terminate at an arbitrary non-singular final state (t^ — > oo). 
No claim of originality is made for this " stochastic interpretation" of the fifth-time 
action. The equivalence of Langevin evolution to the fifth-time action, for ordinary 
bounded actions, was shown by Gozzi in ref. |Tj|. The idea of fixing the final 



state of Langevin evolution for bottomless actions is due to Claudson and Halpern 



13], in connection with Nicolai maps for QCD4, and it was noted by Giveon et. 



al. JTJJ that Langevin evolution with fixed initial and final states is equivalent to 
the fifth-time stabilization prescription of ref. |]]. Section 2 contains an exposition 
of this reasoning, which will then be extended, in section 3, to the special case 
of quantum gravity. It will be shown that the stochastic interpretation fixes the 
integration measure, and suggests a particular choice of operator-ordering, in the 
fifth-time action. 

In section 4 the results of a Monte Carlo simulation of stabilized, latticized, 
Einstein-Cartan gravity are reported. The latticization is similar to that proposed 
by Menotti and Pelissetto in ref. ||15|| , in that local 0(4) symmetry is exact, but 
diffeomorphism invariance is broken by the lattice regularization. The simulation 
is carried out on very tiny lattices: 2 4 x 4 and 3 4 x 4; so obviously the results must 
viewed with caution. Nevertheless, the outcome is interesting. It is found, for a 
range of positive and negative cosmological constants, that the system is always in 
the broken phase < det(e) >^ 0. Moreover, in the neighborhood of A = 0: the 
average curvature approaches 0, the volume per lattice site diverges, and there is 
a peak, possibly divergent, in the negative free energy/lattice site —F. The peak 
in — F near A = and infinite volume suggests that, if the cosmological constant 
becomes dynamical by some mechanism, then its probability distribution is peaked 
near A = 0, as proposed in ref. 0]. Recently Carlini and Martellini [HJ have argued 
that the Coleman mechanism for a vanishing cosmological term is realized in fifth- 
time stabilized quantum gravity; this could be related to the numerical results found 
here. 

An alternative approach to formulating and simulating Euclidean quantum grav- 
ity is the method of summing over simplicial manifolds. In this case the action is 
just the Euclidean Einstein-Hilbert (Regge) action, and it is the lattice regular- 
ization (restriction to simplicial manifolds) which provides a cutoff on the number 
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of manifolds with large positive or negative curvature |T7[. This lattice cutoff on 
large curvatures makes numerical simulation possible, and interesting Monte Carlo 
results have recently been obtained in 4-dimensions [17-20]. It remains to be seen 
whether a theory whose stability is based on the lattice structure is actually inde- 
pendent of the lattice structure, in the sense of universality Also, as noted above, 
the continuation of Minkowski space gravity to Euclidean space is not necessarily 
obtained by replacing the Minkowski space Einstein-Hilbert action by the Euclidean 
Einstein-Hilbert action. Some additional remarks are contained in section 4. 



2 Stochastic Quantization of Bottomless Actions 

This section reviews a selection of ideas found in ref. [12-14], which provide a 
stochastic interpretation for the fifth-time action prescription of ref. [1]. 

What does it mean to quantize a bottomless action S[4>] in Euclidean space? 
Clearly, if the Boltzman factor exp(-S) is non-integrable, the path-integral formu- 
lation cannot be used directly. However, in statistical physics the Boltzman factor 
is just a representation of the probability distribution obtained by random Brow- 
nian motion. From this point of view it is the Brownian motion, as described by 
the Langevin equation, which is fundamental. Quantization based on the Langevin 
equation, known as "stochastic quantization", gives the following prescription for 
thermal averages of operators Q[fa\: 



<Q>= lim — / Dqfx, -T <t 5 < T) Q[<f>(x, 0)}exp[- [ d 5 x if/Ah] (1) 

T^oo Z 5 J J-T 

where 

d 5 <p = -- + v 

0(x,-T) = fa(x) (2) 

In stochastic quantization, the role of the classical action S is to supply the drift 
term in the Langevin equation. For a given noise configuration r], the Langevin 
equation is solved starting from an initial field configuration fa(x) at t 5 = —T. As 
T — > oo, the thermal average should be independent of the initial configuration. 

Applied to the problem of bottomless actions, the Langevin approach does not, 
at first sight, seem much of an improvement over path- integral quantization. Even if 
the initial configuration fa is a local minimum of the action, for large t 5 the field will 
almost always be driven into the bottomless region of the action, evolving towards 
a singular configuration as — > oo. The qualifier in "almost always" is important, 
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however, since there do exist noise configurations 7/(x, £5) for which the field remains 
in the vicinity of a stationary point SS/S(j) = throughout its Langevin evolution. 

The proposal, then, is as follows: Let us retain the prescription of stochastic 
quantization that S[<f)], bottomless or not, supplies the drift term of Langevin evo- 
lution, but average only over noise terms which leave the field in a fixed non-singular 
final state 0/(x) at £5 = T. In other words, 

<Q> = Hm i- /^(x,-T<t 5 <T)5[0(x,T)-0 f (x)]g[0(x,O)] 

xe[- [ T d 5 xr] 2 /Ah] (3) 



where, again, 

d 5 d> = -- + v 

0(x,-T) = &(x) (4) 

If we can show that thermal averages < Q > are independent of the choice of 
both initial <pi and final <pf configurations, as T — > 00, then we have obtained a 
version of stochastic quantization which, for bounded actions, is equivalent to the 
usual version, but which is also meaningful for bottomless actions. I will refer 
to the standard version of eq. (|I]), with a fixed initial state, as "stochastic (I)" 
quantization, and the proposal of eq. (0), with fixed initial and final states, as 
"stochastic (IF)" quantization. 

Following ref. JTII O], let us change variables in the functional integral of eq. 



(0). Writing the Langevin equation as 



8S 

r] = d 5 (f) + — (5) 



we have 



Z 5 = J Dr] 6[<f>(x, T) - <pf(x)]exp[- J* d 5 x r] 2 /Ah} 



0f 5r) 
D<p det\—]exp 

<i 0(f) 

<t>i 5ri 
D(f) det[—]exp 



xe 



-(S[<t>f]-S[<t>i])/2h 



(6) 



Working out the Jacobian, 
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det[ 



<50(x') J 
Gfet[<9 5 5 5 (x - x') + 



5 2 S 



x)^(x') |ts 



cfet[<9 5 (5 5 (:r-a;') + fl(t5-4 



<5 2 S 



det[d 5 ]det[5 5 {x - x') + 9(t 5 - t' 5 



<50(x)<50(x') |t5 
5 2 S 



det[d 5 ]exp 
det[d§]exp 



Trln[5 5 (x-x') + 9(t 5 -t' 5 ) 



^(x)^(x') |t5 J 
<5 2 S 



1 /rf 5 , ^ 
2J 84? 



where #(0) = |, and the meaning of 



5 2 S 



(7) 



W(x')| t5 

is to carry out the D=4 dimensional functional variations, and then replace 0(x) by 
0(x, £ 5 ). Note that only the first term survives in the expansion of the Trace log; all 
other terms in the trace vanish due to the time-ordering enforced 9. Substitution 
of (|7|) and (g) into (]3|) gives 



< Q > = lim £ty(x, -T < £ 5 < T) Q[0(x, 0)] 

or, with a rescaling of £5 t^/2Ti 



(9) 



<Q> = — J D0g[0(x,O)]e" 55 
z 5 J 

s 5 = / A[ I(^ + § L ( |r- 



1 5 2 5, 



(10) 



■<ty' AhS(f> 21 

It is now easy to prove that < Q > is independent of the initial and final states 
and 4>f, since 



<Q> 



lim 



L-mm ^ n 



-{E n +E m )T/h 



En^n[4>f}n[^}e- 2E " T/h 



< ^o\Q\^o > 
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where the E n are the eigenvalues, ^ n the eigenstates (\l/o the ground state), of the 
Hamiltonian corresponding to S% 

f J4 r 1 52 1 ,8S. 2 1 PS, 

K-J^-w + wty-*^ (12) 

which is known as the Fokker-Planck Hamiltonian. Note that whether S is bot- 
tomless or not, the potential term in H 5 is bounded from below, and therefore has 
a well-defined ground state. Since < Q > depends only on the ground state of if 5, 
and not on the initial or final configurations <fii, <pf, we have shown that stochastic 
(IF) quantization, like ordinary stochastic (I) quantization, does not depend on the 
the choice of the initial/final configurations. 

There is one way that the reasoning above could have failed. This would be if 
\&o[0/] = 0, e.g. if 0/ or its derivatives were infinite. Such singular configurations 
would be obtained in Langevin evolution for bottomless actions, as £5 — > 00, if the 
final state were unconstrained. Thus, the only slight restriction we need to make on 
the final configuration is that it is non-singular, ^o[0/] 7^ 0; in fact, this restriction 
is implicit in the choice of initial configuration as well, even for ordinary stochastic 
(I) quantization. 

Equation ( |TUD is the "fifth-time" stabilization prescription proposed in [I], which 
is now seen to be equivalent to stochastic (IF) quantization. It is easy to see from 
([10|) that the H — > limit enforces the classical equations of motion 5S/5<ft = 0, and 
it was shown in that the perturbative expansion of ( |i~0|) reproduces the naive 
perturbation expansion, to all orders in any expansion parameter, generated from 
Taylor expanding e~ s in the usual functional integral.0 

Writing the ground state \l/o of in the form 

^0 = exp[-S eff /2h] (13) 

fifth-time stabilization can be regarded as replaced the bottomless action S by a 
bounded action S e ff, i.e 

<Q> = <^ |Q|*o> 

= l|^(x)Q[0(x)]e^/// 7i (14) 

which has the same classical equations of motion, and the same formal perturbative 
expansion, as S. 

The Yang-Mills Vacuum as Stabilized Chern-Simons Theory 

2 The equality of perturbative expansions holds providing S is stable at zeroth order, which is 
not the case for quantum gravity, c.f. ||]. 
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A curious example of the 5-th time approach is provided by Yang-Mills theory. 
The discussion in this subsection is based on the Nicolai map for QCD4 found in 
ref. [13], but with emphasis laid on the fact that QCD^ itself can be regarded as a 
"fourth" -time action. 

It has long been known [|2TJ that there is an exact solution of the Yang-Mills 
Schrodinger equation in temporal gauge 

J d3 4 [ -w +1 2 F ^ =E * (15) 

with energy E — 0; this is the Chern-Simons state 



*[A] 



exp 



J d 3 xeijk[A^djAl —e a bcA^A b jAl 



exp{cs[A]/2} 



(16) 



But because this state is non-normalizable (the Chern-Simons action is bottomless), 
it must be rejected as a physical state. In fact, due to asymptotic freedom, the true 
Yang- Mills vacuum can be expected to look something like the abelian ground state 
for high frequency fluctuations 



abelian 




[A] = exp 



1 



d 3 xd 3 y Fij 



37r^ j \x — y 

while for low-frequency fluctuations, it was argued in ref. 
the true Yang-Mills vacuum has the form 



■&o[A low ] « exp[-fi J d 3 x TrF 2 ] 



(17) 

see also f22[) that 
(18) 



This expression for the low-frequency vacuum has since been verified by Monte- 
Carlo simulations [23 1 . 

Now although the Chern-Simons state is non-normalizable, it is still an exact 
solution of the Yang-Mills Schrodinger equation, and an interesting question to 
ask is what does the stabilized theory looks like. The answer is quite remarkable: 
Stabilized Chern-Simons theory is the true Yang-Mills vacuum). This is easy to see, 
since the rule is to replace the unbounded distribution exp{cs[y4]} by exp[—S e ff], 
where \1/q — exp[— S e ff/2] is the ground state of the Fokker-Planck Hamiltonian 



d A x 



1 5 2 1 Scs[A] 2 lS 2 cs[A] 



+ 



(19) 



2SAf 8 V 5Af ' 4 5Af 
which in this case turns out to be the Hamiltonian of D=4 Yang-Mills theory in 
temporal gauge, seen in eq. ([15]). As a consequence, stabilized Chern-Simons 
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theory is simply the Yang-Mills vacuum in temporal gauge. Likewise, the "fifth- 
time" action (really a "4-th time" action, since cs[A] is 3-dimensional) is just the 
D=4 Yang-Mills theory. From this fact, it is not hard to see that all (not just 
equal-time) correlators of D=4 Yang-Mills theory can be derived from stochastic 



(IF) quantization of the Chern-Simons action, as shown in ref. |]I3"|| , and in detail 
in ref. |24 . 



3 Stabilized Gravity 

To apply stochastic (IF) quantization to gravity, we must generalize the formalism 
somewhat to allow for field-dependent supermetrics. Denote the fields, which may 
represent the metric, tetrad, or spin-connection, by g N , the supermetric by Gmn, 
and the supervielbien by E^. The corresponding Langevin equation is 3 

d5g M = _ G MNjS_ + E M rl A {2Q) 

where 



< V A (x,t 5 )r 1 B (x?,t' 5 )> = 2h5 AB 5\x-x.')5(t 5 -t' 5 ) 

Gmn = E M E^ (21) 

The analogue of @ becomes 



Z h = J Dg N det[^} exp - J d 5 x [^G MN d 5 g M d 5 g N + ^G 



L Sg 

Xe -!(S[g?]-S[ ff f])/2ft 



Sg M Sg 



(22) 



and from the Langevin equation 



V - E M {d 5 g +G 



Sg 1 



(23) 



we have 



Srj A (x) 
5g L (x') 



E 



M 



5f I d 5 5 5 (x-x') + 



6 



5g L (x') 



G 



MN I 



5S 



5g N (x) lu 5g L (x') 



3 This form corresponds to the (J, = a = case of the Langevin equation for gravity in ref. |2E 
The most general case will not be considered here. 
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SS 



x, r 



, ^ T V(x,r)} 



5g L (x' 

leading to the determinant 



(24) 



det[ 



<iet[95](iet[_E , ]exp 



5s 



x,r 



5^(x) 



Sg M (x, t 5 ) 



(25) 



The derivatives of the supermetric and supervielbein in ( p5|) depend on a choice 
of stochastic calculus. In the Ito calculus, all contractions between G MAr (x, t) and 
r)(x.',t) at equal (fifth) time t are taken to be zero, e.g. 



<^(x,t)^(x,tX(x',t)^(x',t> 
= < <(x,tX(x',t) >< i7 A (x,f)i7V,f) > 



(26) 



This condition can be achieved by defining 



E»(x,t) = E%[g(x,t-St)] (27) 

where St is an infinitesmal fifth-time displacement, so that E% (x, t) is indepen- 
dent of T](pc,t). This fifth-time displacement does not affect invariance under four- 
dimensional, ts-independent, diffeomorphisms. Extending this prescription also to 
the supermetric in (||), i.e. G MAr (x, t) = Ef (x, t)Ef (x, t), results in a determinant 



Sti 

det[— -] = det[d 5 }det[E}exp 

Sg 



- d 5 xG 



MN, 



S 2 S 



x 



2.1 ->Sg M {K)Sg N {*)^ 
Finally, substituting fl28|) into fl22|), we obtain the fifth-time action formulation 



(2* 



-5 5 /7i 



<Q> = i/l>^de«[JS?ig[/ r (x > 0)] C - 

S, = jd*x[\G MN d h g M dsg» + \G* ^ §gN 



ss ss 



Sg M Sg N 



(29) 



We now apply this formulation to the Einstein-Hilbert and Einstein-Cartan 
actions. The Einstein-Hilbert action is 
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Seh = ~\ I d 4 x^g-R (30) 



K 



where n 2 = 167rG. Expanding g^ u = 5^ u + nh^, the action to zeroth-order is 

S° = J^ h^P)p 2 [\p {2) - lP {0 - s) U«M-p) (31) 

where p(° _s ) and are transverse spin-2 and spin-0 projection operators ||26|| . 
In this case the 10 independent fields g A just correspond to the metric components 
Q\iv (/^ > z/ )- The supermetric Gmn is defined implicitly from 



<)<f = I d*xG MN (x)6g M (x)6g N (x) 
' d A x G^{x)Sg^(x)5g al3 (x) 



It is required that the arbitrary constant c in the DeWitt supermetric G^ VOL ^ be 
constrained to c > — |; since otherwise det(G) < and we cannot construct a 
supervielbein. This would break the link between stochastic (IF) quantization and 
the 5-th time action (in fact, S$ would no longer be bounded from below). Applying 
(p9|) to the Einstein-Hilbert action (|30"D, one finds 



S, = - I d b x 



4 

-2hn 2 G 



- 2 G^d h g, v d hgaP + n 2 G-J a ^^- 

! S 2 S 



\ J d 5 x [^G^d 59flu d 5 g aP + l g G-^(R^ - l -g^R) 



4 

x(K*-lg<*0R)-hPy/gR 



(33) 



for the corresponding 5-th time action, where f3 is a singular constant. This ac- 
tion was derived in ||, except that the operator ordering in the last term of S5, 
as well as the functional integration measure, was left undetermined. Stochas- 
tic (IF) quantization determines the integration measure to be the (D=4 dimen- 
sional) DeWitt measure det(E) = \J~G = const, while the "retarded supervielbien" 
Ejtf(x,t) = E^[4>(x,t — 8t)] suggested by the Ito calculus determines the order- 
ing of supermetric and functional derivatives shown above. This ordering will be 
particularly convenient for the Einstein-Cartan theory. 

In ref. [|[ it was shown how to calculate the stabilized 4-dimensional action S e ff 
perturbatively, starting from S5. At linearized level 
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h,Mv 2 \\p {2) + ^ (0 - s) l K,(-p) (34) 

which, like S° above, is transverse; only the sign of the "wrong-sign" conformal 
mode has been flipped. S e ff will be non-local at higher-orders, but it is guaranteed 
to have the same classical equations of motion as Seh- 

The problem with (0), which is a shortcoming of the 5-th time approach in 
general, is that S5 contains higher-derivative terms, in this case proportional to R 2 . 
This means that reflection positivity for reflections across the "ordinary" time (24) 
axis cannot be guaranteed, which is problematic for continuing S e ff to Minkowski 
space.^ The presence of higher derivatives in S5 can be traced to the 2nd derivative 
terms in S. For this reason, it was suggested in || to stabilize the Einstein-Cartan 
theory, which, like the Chern-Simons theory discussed in the last section, contains 
only first order derivatives. The "fifth" -time action corresponding to D=3 Chern- 
Simons theory is D=4 Yang-Mills theory, which is certainly reflection positive across 
any axis. As the Einstein-Cartan theory has a tensor structure similar to that of 
Chern-Simons theory, there is reason to expect that the stabilized version is also 
reflection positive. 

The Einstein-Cartan theory in 4 dimensions has the action 

Sec = / e ^e a A eb A ( dujCd + uCf A ( 35 ) 

In order to write down a single Langevin equation for this system, it is necessary 
to have the tetrad and spin connection u'r in the same multiplet g L , and it 
is convenient to rescale these fields so that they have the same dimensions. The 
dimensional quantities in the theory are k, h, so we rescale 



1 

e = — e 



K 



UJ 



\Fhuo (36) 



with supermetric defined implicitly by [] 



5g 2 = 6e 2 + 6u 2 

= Jd'x^[5e^el + 5Cof5u a u b } (37) 

4 Reflection positivity is not ruled out either, since S5 certainly has this property for bounded 
S despite the higher derivative terms. 

5 This is not the most general possible e — ui supermetric. Only the simplest version of the 
stabilized theory will be considered here. 
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Applying this supermetric to the general formula (Jffi), and noting that in this case 
(as in Chern-Simons theory) the singular term 



G 



MN 



5 2 S 



5g M 5g N 







(3* 



vanishes, the stabilized Einstein-Cartan theory (including a cosmological term 
/ d A x\det(e)) is found to be 



< Q[e,oj] >= — I DeDoo 



Y[det 10 (e) 

x,t 5 



Q[e(x,0),Lu(x,0))e~ S5/h (39) 



where 



K 



is the 5-th time action and 



t;,t; ct (cW w + 2«k 



(40) 



R = du + lo A uo 

T = de + uAe (41) 

are the curvature and torsion two-forms respectively. The fifth-time action ( fTO|) con- 
tains no higher-derivative terms and can be shown, in the lattice version discussed 
in the next section, to be reflection positive across the x^ axis. 

The perturbative expansion of (|4"0"D , at A = 0, is an expansion around the clas- 
sical R = T = solution: 

el = 8 a m + Kb; 
< = ujf(e) + K 2 VhQf (42) 

where 

uf = \e va {d,e h v - d v e\) - \e v \d^ v - d^) - ^e° b (d p e ac - d CT e pc )e< (43) 

is the zero-torsion spin-connection. Then the part of S§ which is zeroth-order in U 
is 
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I / d *Vg 



-g^d^el + 4(-RlR a a g^ -XR + k 2 X 2 ) 



(44) 



where = rfuJ + u; A u7. From (f|3) we see that torsion propagates only at loop level. 
The D=4 stabilized action S e ff was calculated at zeroth-order in k (for A = 0) in 
ref. HI, and it was found that 



d^p 



h^uP 2 



Ip(2) i Ip(O-s) 

4 2 



a/3 



(45) 



where fy^ = + b Ufl is the symmetric part of the tetrad. This result is identical 
to the zeroth order obtained for the stabilized Einstein-Hilbert action. 



4 Numerical Simulation 

There is no known lattice action for general relativity which is exactly invariant 
under a continuous symmetry group analogous to diffeomorphisms. The best one 
can do for the Einstein-Cartan theory is to preserve the invariance under the local 
Lorentz group (0(4) in Euclidean space), and hope that the diffeomorphism invari- 
ance of the continuum action can somehow be recovered at a fixed point. Here 



we follow Menotti and Pelissetto |15|] in introducing a hypercubic lattice with link 
variables 



Ufj,(n) = exp[aP b r b (n) + -aJ bc u b ^(n)] 

= exp[aP b e b ^(n)} exp\^aJ bc u b ^(n)} 

= exp[P b el(n)} exp[^J bc uj b ;(n)} (46) 

where a is the lattice spacing, and from here on we set H — 1. P a and J ab are the 
generators of the Euclidean Poincare group in the four- dimensional spinor repre- 
sentation 



Jab = ^[7a,7b] 

Pa = ^7a(l+75) (47) 



Define the plaquette variable 
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U^(n) = U^(n)U v {n + f i )U- 1 (n + u)U- x (n) (4* 



and lattice curvature, torsion, and tetrad 



KM = ~Tr{J ab [U, v (n)-U^(n)}} 



'"' ' 2 
1 

4 



T;M = --Tr{K a [U, v (n)-U v ,(n)]} 



el = --Tr{KJJ^n)^U?(n)} (49) 

where K a = —^j a (l — 75). Under 0(4) gauge transformations, the lattice tetrad, 
curvature, and torsion transform like vectors and tensors in the latin (local frame) 
indices. 

In terms of these variables, the lattice version of the Einstein-Cartan action 
becomes 



(50) 



which is invariant under local Lorentz 0(4) gauge transformations, but still bot- 
tomless. The lattice cosmological constant A is expressed in units of the Planck 
length, i.e. A = K A \ contmuum ^ note that the lattice spacing has dropped out of the 
action. In terms of the same variables, the latticized version of the Einstein-Cartan 
5-th time action (EOT) is 



1 



S 5 = jJ2\det(e)\ -eM<(n){T^(n)T v %(n)+R%(n)Rl b 5 (n)} 

+HKKAe p A< - a<« + a 2 } 



"pv "pa^d^d^b^c '*^"fj,u 
" fiv pa \ 



+eT° u T b ($ ab e^ c + 2e£e p b )e v d e° d (51) 



where 



(52) 



_ _05_ 

and 05 is the lattice spacing in t 5 direction. The lattice vectors n now label sites on 
a D=5 dimensional lattice, and is inverse to the matrix e a . Defining U^{n) = 1, 
the curvature and torsion R^ 5 (n) and T® 5 (n) are also obtained from fl49|). The 
quantities in fl49|) are dimensionless; the scalings which give the usual curvature, 
torsion, and tetrad in the continuum limit are 



Tjab(continuum) f? ab 

V> v ~ a 2 V- u 
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rpab(continuum) rpab 

1 nu — a 2 H» 



a(continuum) c a 

e M - -e^ [06) 



The action ([y]), together with the integration measure 

/ n^t 10 (e)n^H (54) 

J n n 

was used for the numerical simulation. Note that the lattice spacing a in the 
1 — 4 directions scales out of the action completely, and all quantities which, in the 
continuum, are diffeomorphism invariant, appear on the lattice in units of k (c.f. 
13), e.g. 



Itfxy/gR = K 2 Y.\det{e)\R a y a el 

J n 

I d 4 x-s/g = k a \det{e)\ 

J n 



5s 2 n>n+s = K 2 <ey v >5»5» (55) 

The expectation values of such quantities are therefore expressed, on the lattice, in 
units of the Planck length k rather than lattice spacing a, and depend only on the 
parameters e = a<j/2K 2 , and A. The fact that lattice spacing a drops out of both the 
action and the expectation values is a remnant of diffeomorphism invariance in the 
continuum, and has led to the speculation that lattice quantum gravity, like string 
theory, has a kind of built-in minimum length 5s 2 ~ k 2 . 

The final step in the lattice formulation above would be to symmetrize fl5"T| ) 
with respect to all 7r/2 rotations around the axes; the resulting symmetrized action 
can be shown to be reflection positive fl5| . In the interest of minimizing computer 



time however, only the unsymmetrized action ([H]) was used for the Monte Carlo 
simulation. 

Computer simulation of latticized Einstein-Cartan theory is rather lengthy, even 
for very small lattices. The link variables are 4x4 matrices on a 5-dimensional 
lattice, and because of the plaquette-plaquette structure of the action, a single link 
update requires on the order of 10 4 floating-point multiplications. For this reason 
the simulation has only been carried out on tiny lattices of dimensions 2 4 x 4 and 
3 4 x 4 (4 spacings in the £5 direction), which required, despite the small lattice sizes, 
a total expenditure of 100 Cray YMP hours. 

There is a two-parameter (e — A) space of couplings for 65. One cannot vary e 
at fixed A = because the lattice does not seem to thermalize at A = 0, for reasons 



6 Disappearance of the lattice spacing a in favor of k is found also in the Regge lattice 
formulation. 
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discussed below. Instead I have fixed e somewhat arbitrarily at e = ~, and carried 
out Monte Carlo simulations for various values of A. 
Figures 1-3 plot the average curvature (Fig. 1), 



< R > 



< £ \det(e)\R > 
<J2\det(e)\ > 



(56) 



the volume per lattice site (Fig. 2), 



<Va>=^- T — <El det ( e )l> 



(57) 



sites 



and the derivative of the free energy per lattice site (Fig. 3) 



dF 
~d\ 



sites 



1 



<Y,\det(e)\(-.5R + A) > 



n 



(58) 



versus the cosmological constant A in a range —5 < A < 20. Squares denote data on 
2 4 x 4 lattices, at A = —1, .25, .5, 1, 2, 5, 10, 20, crosses denote data on 3 4 x 4 lattices 
at A = —5, —2, —1, .5, 1, 2, 5. The longest runs were at the smallest values of A, e.g. 
300 thermalizations and 1600 data-taking iterations at A = 0.5 on a 3 4 x 4 lattice. 

From Fig. 1, we find that < R > depends linearly on A, with < R >^ 1.6A, 
as compared to the classical value R = 2A. The modest discrepancy in slope could 
have various origins: renormalization of A, the smallness of the lattice, and the 
effect of the hypertoroidal lattice topology. 

Figure 2 shows a clear divergence in volume/site at A = 0. Runs at the value 
A = did not converge to a finite value of |det(e)|, instead showing a steady increase 
in volume with number of iterations. Apart from the divergence at A = 0, there 
are two other noteworthy facts concerning < det(e) >: First, it is found that for all 
values of A < 10, < det(e) >=< |<iet(e)| >, while for the largest value of A = 20, 

< det(e) > was only 4% less than < |det(e)| >. From this we conclude that in 
the range of A studied, the system is in the broken symmetry phase expected for 
quantum gravity, and it appears that fluctuations which change the sign of det(e) 
are rarely generated.f] The second point is that the volume/site decreases smoothly 
away from A = for both positive and negative A. This fact, and the fact that 

< R >oc A, show very clearly the effect of stabilization, since the cosmological term 
in the Einstein-Cartan/Hilbert actions is bottomless for A < 0. In the fifth-time 
Einstein-Cartan, changing the sign of the cosmological constant (leaving Newton's 
constant fixed) is equivalent to changing the sign of Newton's constant (leaving the 

7 Strictly speaking, transformations such as inversions — > — x 4 , U± — > U^ 1 , which take 
det(e) — > —det(e), are invariances of the symmetrized, reflection positive lattice theory. For 
the unsymmetrized action of eq. (|5l|), used for the numerical simulation, invariance under such 
transformations is only approximate. 
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cosmological constant fixed). In the Regge and simplicial manifold approaches, the 
behavior of average volume and curvature changes drastically upon changing the 
sign of Newton's constant [17-20]. This is not the case in the stabilized version of 
the Einstein- Cart an theory, where the metric fluctuates close to zero curvature for 
small A of either sign. 

In connection with simplicial manifolds, it has been noted that the lattice cut- 
off stabilizes the Euclidean Einstein-Hilbert action, since the entropy of simplicial 



manifolds with large curvatures is small [0. If continuum quantum gravity would 
then be stable due to the measure, the 5th-time procedure cannot change anything, 
and is simply a rewriting of the original 4-dimensional theory. (This is because if 
exp(-S) is normalizable, then it is necessarily the ground state (squared) of H 5 ; 
this means that S e ff = S, where S is the Einstein-Hilbert action.) On the other 
hand, it may be that any stability induced by the lattice cutoff in the D=4 theory 
is spurious or non- universal, in which case the 5th-time action is a much better 
starting point for both numerical and analytical work. Perturbative calculations 
in the stabilized Einstein-Hilbert theory do indicate, in fact, that S e ff ^ S, and 
indeed that S e ff is non-local ||, which is in accord with the findings of ref. ||, |9|. 

Finally, we see in Fig. 3 that the slope of the free energy is large, and possibly 
singular, in the neighborhood of A = 0. Although the data is nowhere near sufficient 
to allow a numerical integration, it does appear that the negative free energy (and 
therefore Z 5 ) must have a peak, and might even be divergent, at A = 0. This peak 
is intriguing, because if the cosmological constant becomes dynamical in some way, 
as suggested by Baum, Hawking and Coleman [f|, then Z 5 (X) has the interpetation 
of a probability density for the cosmological constant, i.e. 



P(A) = N j Dg e ~ 5e « [9 ' A] 

= jDgil>*\g,\]i/>\g,\)<xZ 6 (\) (59) 

A sharp peak in the (negative) free energy at A = is therefore an explanation of 
the smallness of the cosmological constant. It is too much to claim that the data 
shown here, obtained on tiny lattices from an action which breaks diffeomorphism 
invariance, is strong evidence for the vanishing of the cosmological term in quantum 
gravity, but at least the data does seem to support this idea. Actually the evidence 
for a sharp peak in — F is better for A — > + than for A — > 0~. Conceivably, since 
the volume is divergent at A = 0, Z§ itself may be discontinuous at this point, with 
the probability peaked at A = + . As usual, one would like to have more data, on 
larger lattices, in the region of interest. 

I have made no attempt to search for an ultraviolet fixed point in the e — A 
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parameter space; all computations at various A were made at the constant (and 
rather arbitrary) value of e = \. A fixed point where diffeomorphism invariance is 
restored (if such a point exists) would be signaled by a peak, growing sharper with 
lattice size, in the correllation 

c = <(EVgR) 2 >-<EVgR> 2 () 
<EVa> 

No such peak was observed for e = \ on the 3 4 x 4 lattice; the value of C xs 11.5 
was almost constant, varying by less than 10%, over the full range of A. A very 
interesting question is whether a peak in correllation C develops as e is varied; this 
is an issue where much additional numerical work is required. 
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Figure Captions 



Fig. 1 Average curvature vs. cosmological constant, in units of the Planck length. 
Squares denote data on 2 4 x 4 lattices, crosses denote data on 3 4 x 4 lattices. 

Fig. 2 < ^fg > (volume/site) vs. cosmological constant. 

Fig. 3 Derivative of the free energy/lattice site dF/dX vs. cosmological constant. 
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